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Abstract 

Several determinants with gamma functions as elements are evaluated. These kinds of determinants 
are encountered, for example, in the computation of the probability density of the determinant of 
random matrices. The s-shifted factorial is defined as a generalization for non-negative integers of 
the power function, the rising factorial (or Pochammer's symbol) and the falling factorial. It is a 
special case of polynomial sequence of the binomial type studied in combinatorics theory. In terms 
of the gamma function, an extension is defined for negative integers and even complex values. 
Properties, mainly composition laws and binomial formulae, are given. They are used to evaluate 
families of generalized Vandermonde determinants with s-shifted factorials as elements, instead of 
power functions. 

1. Introduction 

This work has been motivated by studies of the probability density of the determinant (PDD) of random 
matrices [1-3]. The method used, and sketched in section 5, is to compute the Mellin transform of the 
PDD. In many cases it turned out to be a determinant with gamma functions as elements. One aim of this 
work is to evaluate some of these determinants and more generally determinants with shifted factorials (or 
Pochhammer's symbols) as elements. 

We define in section 2 the s-shifted factorial (z)s-n, equation (2.1), as a generalization for non-negative 
values of n of the power function z", the rising factorial {z)n, equation (2.3), and the falling factorial [z]„, 
equation (2.4); both the names and the notations of these last two objects are not well established, see [4- 
10] ^ . As a function of z, the s-shifted factorial is a special case of the polynomial sequences of the binomial 
type studied mainly in the calculus of finite differences and combinatorics; see in particular [6, 7] and for a 
wide bibliography [8, 9]. Expressed in terms of gamma functions, the s-shifted factorial can be extended to 
negative values and even complex values of n. The s-shifted factorial provides a compact formulation which 
emphasizes similarities and connections which exist between the power function and the shifted factorials: 
multiplication laws, Pascal triangle property, generating function and binomial formulae. 

It is shown in section 3 that Vandermonde determinant with {zj)s;i instead of {zjY as elements is still 
equal to the usual Vandermonde determinant. Other determinants with the inverse of a s-shifted factorial, 
or the ratio of two s-shifted factorials, as elements are also evaluated, both for positive and negative values 
of the index i. Using the relations between the s-shifted factorial and the gamma function, or the binomial 
coefficient, to each determinant evaluated in section 3, it corresponds a determinant in section 4 with elements 
expressed in terms of gamma functions. Finally, some applications of these determinants are given in section 
5: evaluation of the PDD of random matrices and also, a possible application to Stieltjes moment problems 

^ See [4] Pochhammer's symbol (z)„ := 2(2 + 1) • • • (z + n - 1) 6.1.22, [5] (z)„ -.^ z{z + 1) ■ ■ ■ {z + n - 1) 
p. xLiii, [6] ^ factorielle z descendante d'ordre n' (z)„ := z{z — 1) • • • (z — n + 1) [4f], ^ factorielle z montante 
d'ordre n' or Pochhammer's symbol n < z >n'-— z{z + 1) • • • (z -|- n — 1) [4g], [7] lower factorial (z)„ := z(z — 
1) • ■ • (z — n-|-l) (1.1), upper factorial z'"' := z(z-|-l) • • • (z+n — 1) (1.2), [8] section 5 falling factorial sequence 
(z)„ := z(z — 1) • ■ • (z — n + 1) 2.1, rising factorial sequence < z >„:= z(z + 1) ■ • • (z -|- n — 1) 3.1, [9] falling 
factorial of length n [z]„ :— z(z— l)---(z — n+1) 3.2, rising factorial of length n [z]" := z(z+ 1) • ■ • (z-|-n — 1) 
3.4 and III. 2. A and [10] z— := z(z — 1) • ■ • (z — n -I- 1) nth falling power of z and z" := z(z -f 1) • ■ • (z -f n — 1) 
nth rising power of z 3.4.2. 
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arising in connection with the boson normal ordering problciin. As anotlier example of application of the 
binomial formula, the finite sum of s-shifted factorials of arithmetic progression to n terms is evaluated in 
appendix A. Some basic properties of the product of differences and of the Vandermonde determinant are 
recalled, respectively, in appendices B and C. Finally, appendix D illustrates another way to handle s-shifted 
factorials. 



2. Definition and some properties of the s-shifted factorial 



2.1. Definitions and relations between shifted factorials 

With n a non-negative integer, z and s (the shift) some complex numbers, let us define the s-shifted factorial 
by 

n = 

z{z + s) ■ ■ ■ [z + {n — l)s) n = 1, 2, 

For s = 0, 1 and —1, this definition coincides, respectively, with the power function, the rising factorial (or 
Pochhammer's symbol, mainly in hypergeometric theory) and the falling factorial, namely for n nonzero, 



(^)«;n:={^^^^ , ,^ , (2.1) 



{z)o-n = (2.2) 

(^)l;n = (2)n ■.= z{z^\)---{z^n-\) (2.3) 

(^)_i;„ = \z\n ■.= z{z-\)---{z-n + \) (2.4) 

and when n = all these quantities take the value 1. Thereby, the s-shifted factorial allows compact 
expressions which emphasize the similarities between the power function and the shifted factorials. 
For any non- negative integer n, one has 

(2),;„ = (-l)"(-Z)_,;„ (2.5) 
= (^+(n-l)^)_,;„- (2.6) 

For s nonzero, the s-shifted factorials are related to the rising factorial, equation (2.3), by 

z 



(^).;n = s" - . (2.7) 

\S / n 

As a function of z, {z)s;n is a monic polynomial (i.e. the coefficient of the highest power is one) in z of 
degree n, 

{.z)s,n = z^ + s ^"-1 + • • • + (n - 1)! s"-i z (2.8) 

with 0,— s,...,— (n— l)sas zeros. Consequences of these properties in terms of Vandermonde determinants 

are developed in section 3. The sets of polynomials {(z).,;,, , n = 0, 1, . . .} are special cases of the polynomial 
sequences {Y>n{z), n = 0, 1, . . .}, Y>n{z) being exactly of degree n. We are going to use these sequences in the 
way it is done in combinatorics [7-9] ^ . Any polynomial sequence is a basis of the vector space V over the 
complex field of complex polynomials in the variable z. Then, to any two polynomial sequences {p„(z)} and 
{q„(z)} there exist uniquely determined connecting coefficients such that qn{z) = Yllt=Q^n,kVk{z). These 
important coefficients have been widely studied, e.g., 

n n n 

Wn = E^("'^)^'' ^" = ^^(n,fc)Mfc {z)n = J2L{n,k)[z]k (2.9) 

fe=0 fe=0 k=0 

where s{n, k), S{n, k) and L{n, k) — {^'2^n\ / k\ are, respectively, the Stirling numbers of the first and second 
kind [4, 6, 7, 9, 10] ^ and the signless Lah numbers [6, 7, 9, 10] ^ (other relations between z", (z)„ and [z]„ 

^ See, e.g., [7] section 1, [8] section 3, [9] section III.2. 

3 See, e.g., [4] 24.1.3,4, [6] chapter V [5e,f], [7] (1.11-13), [9] 3.24,25 or [10] 2.5.2. 
^ See, e.g., [6] chapter III p. 165 , [7] (1.14) and section 9, [9] 3.24,25, [10] 3.1.8. 
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immediately follow from equations (2.5) and (2.6)). We will see in subsection 2.6 that {(z)s;„,n = 0, 1, . . .} 
has in addition the important property to be a polynomial sequence of the binomial type. 

2.2. Special values 

With k some non-negative integer, one gets 

{0 k = 0,...,n-l 

k = n,n+l,... (2.10) 
[k — ny. 

r fc = 

(fc)„ = (-l)"[-fcU= (fc + n-l)! (2.11) 

Then, for s nonzero, values of {ks)s;n follow from equation (2.7), in particular (s)s;„ = n!s". 

2.3. Relations with the gamma function and definition of the generalized s-shifted factorial 
One has [4] ^ 

(^)-^^^T^^ <-"•'•" (7) 



T{z-n+l) {z-n)\ 



For s nonzero, relations with (2;)s;„ follow from equation (2.7). 

Actually, the relations above can be taken as the definition of {z)s;n in terms of the gamma function. 
Thereby, one extends the s-shifted factorial to negative values, and even to complex values t of n, defining 
the generalized s-shifted factorial by 

r(2. + t) 

W.;*:=«*4^ (2.14) 

with a cut, say, along the negative real axis of the complex s plane, with — tt < args < tt, to ensure a 
single- valued dependence on s, and we choose the determination such that s* = 1 for s = 1. As s goes to 
zero, say, along the real axis, using the Stirling formula [4] ^, r{z) ~ e~^z^~^ {21^)^(1 + 0{l/z)) as — > oo 
in |argz| < tt, one recovers the function 2;*. Furthermore, the definitions of the rising and falling factorials 
are extended by 

{z)t := {z)i.,t [z]t := (^)-i;t . (2.15) 
Then, it immediately follows from the definition (2.14), 

{Z)s;0 = 1 (2.16) 

i^Ut = {-zU;t (2.17) 

by the recurrence formula r(^ + 1) = zr{z), 

iz)s;l = Z (2.18) 

and by the reflection formula [4] T{z) r(l — z) — n/ sin(7r2;), 

, ^ s* sin(7r-) , , ^ , . 

(^)«;* = , , .^^ i' + it- • (2-19) 



i-sY sin(7r(f +t)) 

5.1.21 and 6.1.2 

^ See, e.g., [4] 6.1.17. 



5 See, e.g., [4] 6.1.5, 6.1.21 and 6.1.22. 
<5 See, e.g., [4] 6.1.37. 



Thus for any integer q, 

{Z),,, = {-IY{-Z).s;, (2.20) 

= {z + {q-l)s)_^,^ (2.21) 
which generahze equations (2.5) and (2.6) for any integer, even negative. 
2.4- Multiplication laws 

When the power function fulfils z* z^ = with r and t some complex numbers, it follows from the 

definition (2.14), 

{Z)s,t{z^ts)s,r = {z)s;t+r (2.22) 

and in particular, by equation (2.16), setting r = —t yields 

(^)«;* = TTrh • (2-23) 

[Z + tS)s--t 

This relation generalizes z^ = l/z~* for s = and provides the relation between the s-shifted factorials for 
any integer q and —q, by equation (2.21), 

(z - qs)s;q {z - qs) {z - {q - l)s) ■ ■ ■ {z - s) {z - s)_s;q ' ''^'^^^ 
In terms of binomial coefficients, the multiplication law (2.22) reads 

:)C;X"r)U.) - (:)H-W'(r-:)- 

For a proper choice of determination the power function fulfils w^z* = {wzY- For the s-shifted factorial, 

one has ^ 

{wz),.t = j^{z)^.t (2.26) 

and thus for any integer q, 

iwz),., = uf'{z)^,q. (2.27) 

For w = —1, relation (2.26) corresponds to equation (2.17). For w = k and q = n some non-negative integers, 
iterating the multiplication law (2.22) and from equation (2.27), 

k—l n—1 k—1 Q 

{kz)s-kn = fc'" n IK^ + + i) = IK-^ + ')s;n (2-28) 

where the last equality corresponds to a rearrangement of the factors, both nl + j and £ + jk taking once 
all the kn values 0, 1, . . . , A;n — 1. The equation above can also be obtained from the definition (2.14) and 
the Gauss multiplication formula [4] ^ , r{kz) = (27r)i(i-*=)A;'=^-3 H^Iq r(z + £/fc). Note that, based on the 
reflection formula and the Gauss' multiplication formula, 2sin(7rA;^;) follows the known multiplication law 
similar to equation (2.28) [5] ^ , 

2 sin(7rA;z) = 2 sin('7r(z + -)) . (2.29) 



See, e.g., [4] 6.1.20. 
See, e.g., [5] 1.392 (1.). 
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For a proper choice of determination the power function fulfils (z^Y = z*''. No equivalent general 
relation exists for the s-shifted factorial. Although {z~^)p has no simple relation with let us 

point out the following expression for any integers n > p > 0, by equations (2.24), (2.22) and (2.6), with 
z^O,-s,...,-{n-l)s, 

' -(.+A-,= (-?:^— = ^^^^\";^^^-— (2.30) 



recovering for s = the relation (z^)^^ — z^p — z"^p (z")~^. 

2.5. Generalized Pascal triangle property and s- difference operator 
The multiplication law (2.22) and equation (2.18) yield 

{z)s-t -{z- s)s.t = ts {z)s;t-i (2.31) 
which generalizes the Pascal triangle property for binomial coefRcients, by equation (2.12), 



z + l\ z\ z 
n I \nl \n — l 



(2.32) 



Let us define the s- difference operator Ag on functions f of 2; by 

A J(z) := i{z + s) - i{z) (2.33) 

(this operator must not be confused with the product of differences A„ (z) introduced latter in section 3 and 
defined by equation (B.l)). It follows immediately from equation (2.31), 

K{z)s-t = ts{z + s)s-t-\ A_s{z)s;t = -ts {z)s;t-i (2.34) 

and iterating these formulae, e.g., the first one 

APiz)s:t = [t]psP{z+ps)s;t-p (2.35) 

recovering for s = the expression of ■^z"'. 

2. 6. Generating function and binomial formulae 

With X some complex variable, let Gs-z{x) be the generating function of the s-shifted factorials {z)s-n-, 



GsAx) ■■= Y.{z)s;n ^ , |sa;| < 1 (2.36) 



n=0 



and, using equation (2.7), 

G.;.(a;) = Gi;^(sar). (2.37) 

Now, the generating function of the rising factorials can be obtained directly from the binomial series with 
equation (2.12), 

(1 - x)-^ = ^(-1)" "J X" = ^(z)„ ^ = Gi;,(x) , |x| < 1 . (2.38) 

n=0 ^ ^ n=0 

Therefore, 

Gs;,{x) = {l-sx)^ (2.39) 
recovering for s = the expression Go,^(a;) := X^^g ^ ~ 
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Since the generating function Gs;z{x) above reads as an exponential function F(a;)'^ of z, it satisfies the 
multiplication law [11] 

Ge;z {x) G,;^ {x) = G,;,+^ {x) . (2.40) 
Expanding both sides of this last equation as a power series in x yields, 

{Z + W),.ri = X! ( t ) Ms;n-k (2-41) 

fe=0 ^ ^ 

namely, the s-shiftcd factorial satisfies the binomial formula. The polynomial sequence {(z)s;„,n = 0, 1, . . .} 
which satisfies (z)s;o = 1 and the binomial formula above is said to be of binomial type [6-10] ^'^ . This 
property is shared by many other binomial sequences {p„(z). n = 0. 1, . . .} which have been studied mainly 
in combinatorics using generating function methods and above all efficient operator methods. 

The binomial sequences can be characterized by a generating function which depends exponentially on 
z [9] 11 

n=o 

then pi^\z) is a monic polynomial of degree n in z, the coefficients of which are known as Bell polynomials 
[6] 1^ (indeed, expanding the exponential series above, the term in z'^ reads 2"a;"(l + 0{x))/n\). In the 

case wc consider, g{x) := —s^^ ln(l — sx) and pn^\z) — {z)s;n', for s = 0, g{x) = x and pn^\z) = z". The 
binomial sequences can as well be characterized by the fact [7-9] that the basis operator of the sequence, 
i.e. the linear operator D of the vector space V (already considered in subsection 2.1) into itself defined by 
Dpo(z) := and _Dp„(z) := npn-i(2) for n > 1, is a delta operator, i.e. it is shift invariant, DEa = EaD 
for all complex number o, where Ea is the translation operator defined by Eat{z) := f(a + z) and moreover 
Dz = c ^ 0. In our case, from equation (2.34), Dg := — A_s/s = {l—E-s)/s, where I is the identity operator, 
is clearly shift invariant and DgZ — 1. Indeed, the binomial formula (2.41) for the s-shiftcd factorial can also 
be demonstrated by recurrence from the Pascal triangle properties (2.31) and (2.32). It is true for n = 
and 1. Let us assume it to be true for n, then, 

{z + w)s;n+i = X] ( ) ^^^^''^ {z + ks + w + {n-k)s) 

k=0 ^ ^ 

^ E ( _ 1 ) ^^^^''^ i'^)s\n+l-k + ^{,] {w)s;n+l-k 



k=0 ^ ' fe=0 

n + 1 



By equation (2.12), in terms of binomial coefficients, the binomial formula (2.41) reads 



n 



k=0 



As for the power function, the binomial formula for the s-shifted factorial can be directly extended to p > 2 
variables using multinomial coefficients, 

(E E ;^^T^(-i).;..---(-p).;n.. (2.45) 

7=1 "1 np=o y 

ni + --- + np = n 

1° See, e.g., [6] [6a] and [13c,d], [7] (1.6), [8] section 5, [9] section III.2, [10] 3.4.2 (3.). 

" See, e.g., [9] 3.59. 

12 See, e.g., [6] section III.3. 

1^ See, e.g., [7] section 3 theorem 1, [8] section 7, [9] 3.45. 
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Prom equation (2.5) the following corollary is immediately obtained: 

{Z - W)s;n = E(-l)"-' (l) {Z)s;k M-s;n-k ■ (2.46) 
fe=0 ^ ' 



Although, as already noted, {z ^ {{z)i) ^, the binomial formula can be extended to the inverse of 
s-shifted factorials. Indeed, by equations (2.30), (2.41) and (2.6), 

V [""^ ^ = V [""^ (^ + ("-^»-.;n-fc (^+(^-l)g)_,.fc 

- Ill + 2(n — 1 I ( z + in + (n ~ 1 )s] 

(2.47) 



(z + w + 2(n-l)s) . (z + w + {n~l)s) , 

{z)s;n {w)s;n {z)s;n {w)s;n 

corresponding for s = to 



^ \k J z'' w''-'' z'^w'^ \z w I ^ ' 



k=0 



where, once again, the last equality above does not hold for s nonzero. Similarly, one also gets from equations 
(2.30) and (2.46) 

,.n^k{Aiz)s;k ^t,.n-kfA.. (^+(n-l)5)^ 

{z-w-{n-l)s) 



iw)s;r. 

corresponding for s = to 



(2.49) 



k=0 



Another kind of useful relations is as follows. From equation (2.25), the multiplication law (2.22) and 
the binomial formula (2.41), 

Yl (fc) ['^JP M^;n-k = [n]p X! (fc _p) (^)^;* {w)s;n-k 

= [n]p {z)s-p ( ? _ ^ ) {z+ Ps)s;k-p {w)s;n-k 

k=p ^'^ 

= [n]p {z)s-p {Z + W+ ps)s;n-p ■ (2-51) 

Several of these binomial formulae are used in the next section to evaluate some determinants with 
s-shiftcd factorials as c;lc;nic!nts. As another example of application, the finite sum of s-shifted factorials of 
arithmetic progression to n terms is evaluated in appendix A. 



3. Generalized Vandermonde determinant with s-shifted factorials as elements 



In what follows, n is a positive integer and z is, either a set of complex numbers, or a complex function, 

^■={zj,j = 0,...,n-l} or j z{j) := Zj j = 0, . . . , n - 1 . (3.1) 

7 



Some basic properties of the product of differences A„(z) := no<i<j<n-i(-^i ~ ■^«)' equation (B.l), and of 
the Vandermonde determinant det [(-s^j)*] j j_q are recalled, respectively, in appendices B and C. 

3.1. Expressions for s-shifted factorial with a non-negative index 

Lemma 1. With n a positive integer and s some complex number, the generalized Vandermonde determinant 
of s-shifted factorials, still is the product of differences, 

det[(^j)^;»]ij=o,...,„-i=^"(^) (3-2) 
thus, it does not depends on s. More generally, 

det[ni(^,)], ^^_^=AA„(z) (3.3) 
where Ili{z) are n linearly independent polynomials in {z)s;. each of degree less than n and defined as follows: 

n-l 
fe=0 

In particular, with hi some complex numbers, one has 

det [{bi + zj)s;i\ ij^o,...,n-i = ^"(^) • (3-5) 
Finally, with t some complex number, 

n-l 

det [{Z3)s;t+i]ij^0,...,n-l = {I[^^j)s;t} A„(z) . (3.6) 

Two proofs are given. Based on the properties of the s-shifted factorial, proof 1 expresses the determi- 
nants considered in terms of Vandermonde determinants. Illustrating again the similarities between (z).,;/ 
and z*, proof 2 follows the same steps as a usual way of computing the Vandermonde determinant, equation 
(C.l). 

Proof 1. The s-shifted factorial (z)s;i is a monic polynomial of degree i in z, sec equation (2.8). Hence, 
formula (3.2) follows from equation (C.3). Note that equation (3.2) still holds for the monic polynomials 
obtained from any generating function defined by equation (2.42). Formulae (3.3) and (3.5) can be directly 
obtained either from equation (C.3) in terms of usual polynomials (e.g., {bi + z)s;i is also a monic polynomial 
of degree i in z) or starting from formula (3.2), by the same arguments as for equation (C.3), in terms of 
polynomials in s-shiftcd factorials (e.g., by the binomial formula (2.41), {bi + z)s-i = J2k=o {k){bi)s;k{z)s;n-k, 
i.e. a monic polynomial of degree i in {z)s;.). Finally, equation (3.6) follows from the multiplication law 
(2.22) and formula (3.5). 

Proof 2. Let Mij := {zj)s;i- The determinant dct[Mi,j]i j=o,...,n-i is not changed if one replaces the row 
TZi by the linear combination TZi — (Mifi/Mi-i^) TZi-i, successively for i = n — 1, n — 2, . . . , 1. Then, by the 
multiplication law (2.22), for i = 1, . . . , n — 1 and j = 0, . . . ,n — 1, 

Mij {zj)s-i - {zo + {i - l)s) {zj)s;i-i = {zj - zo) {zj)s-i^i . (3.7) 

This operation replaces the column Co by zeros except for the row TZq left unchanged. Expanding the 
determinant with respect to Co and taking out the factors depending only on j yield the recurrence formula 
on n, 

D^-nizo, Zn-i) := det[Mij]..^^^^^_^ 

n-l 

= - Zo)^ Ds-n-lizi, ■ ■ ■ ,Zn-l) ■ (3.8) 
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An iteration of this last equation, down to Ds-i{zn-i) = 1, completes the proof. 

The recurrence procedure above makes step by step the matrix (Af^ j)^ j^g^ „_i triangular. Let us 
denote by a superscript the rank of the step in this procedure. At the first step the row i = is unchanged 
while for i = 1, . . . , n — 1, 



7^(1) = 7e Mi^n- 1 

At the second step the rows i = 0, 1 are unchanged, while for z = 2, . . . , n — 1, 



(3.9) 



(1) 



-7?(l) 



i,0 



M, 



(1) 



(1) 



-1,1 



Generically, the final expression of the row i is given by 'R.^^ \ It happens that in the special case Zj := b + js 
with b some complex number and s nonzero, these expressions read 



(1) M, 



Hi. 



i-2,0 



(3.10) 



7^, 



(1) 



(fe + (i-i)s)7ei-i 



Tef ^ = - 2(6 + (i - l)s) 7ei_i + (6 + (i - l)s) (6 + (i - 2)s) 7li_2 



(3.11) 
(3.12) 



/s=0 



(3.13) 



where the last formula can be checked as follows. With Mi^j := [b + js)s;i, by the binomial formula (2.46), 
equations (2.7) and (2.6), 



(3.14) 



which vanishes for i > j, see equation (2.10). Another proof of this identity is given in appendix D.l. Thus, 
as expected, the resulting matrix is triangular and its determinant is the product of its diagonal elements 
= j\. Then, by equation (B.5), 



det [(6 + Js).;.],,,^o,...,n-i = s"^"-'^/' n = ^"(-^^ -^b + js) 
completing the proof of equation (3.2) in the special case zj := b + js. 

Lemma 2. With n a positive integer, s some complex mmiber and Zj ^ 0, — s, . . . , — (n — 2)s, 



det 



(_l)"(«-l)/2 
rn-l/ 



A„(z). 



i)j=o,.--)"-i 11^=0 (•^j)s;"-i 
This formula generalizes equation (C.5). 

Proof 1. With n — 1 > i > and Zj ^ 0, — s, . . . , — (n — 2)s, by equations (2.30) one gets 

i,j=0,...,n-l 117=0 {.^j)s;n-l 



det 



1 



(•2j)s;i 



(3.15) 



(3.16) 



(3.17) 
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Then, changing i into n — 1 — i (i.e. rearranging the rows) on the right-hand side determinant above, lemma 
2 follows from equation (3.5). Note that when s = 0, equation (C.5) for the power function can also be 
derived as above from = z"'~^~^ / z"'~^ . 

Proof 2. This proof of lemma 2 follows the same steps as proof 2 of lemma 1. With the linear combination 
of rows Tli - (Mi,o/Mj_i,o) T^i-i, 



M,. 



1 



1 



''^ ' {Zj)s;i {zo + l)s) Zj {Zj + s)s-i-\ ' 

Then, with Zj ^ 0, — s, . . . , — (n — 2)s, the recurrence formula on n reads 

Ds-n{zo, • . • , Zn-i) := det [M^ , ,,-^o_ „_i 



, n — 1 . 



{zQ)s;n-l n?=l Zj 



Ds-n-\{z\ + S, . . . , Zn-1 + s) . 



(3.18) 



(3.19) 



Iteration of this last equation, down to Dg-i (zn-i + {n — l)s) = 1, ends the proof. 

As in proof 2 of lemma 1, in the special case zj := b + js with b some complex number and s nonzero, 
the determinant can be made triangular in one step, replacing TZi by the linear combination 



Vfc; (-5-2(z-l)s) 



■ 7^fe, 



(3.20) 



s:i—k 



Indeed, with Mij := l/{b + js)s-i and b^O, —s, . . . , — (2n — 3)s, it follows from the binomial formula (2.47) 
and equations (2.5)-(2.7) 



M, 



E 



1 



ykj {-b - 2{i - l)s)^.._^ (6 + js),,k {b+{i- l)s)^.. (6 + js)s;i 



(3.21) 



which vanishes for i > j. Another proof of this identity is given in appendix D.2. Then, the determinant is 
the product of its diagonal elements. 



det 



1 



{b + js)s-i 



n-l 



i,j=0,-..,n- 



-1 ^ ' j-}jb+{j-l)s)^..{b + jsU, 



and, using the multiplication law (2.22), it can be shown by recurrence that for all s 

n—l n—1 

l[{b+{j - 1)S)^.. {b + js)s;j = l[{b + js) s-n-1 

corresponding to a rearrangement of the factors. Finally, by equation (B.5), 



det 



{b + js)s;i 



(_l)"(n-l)/2 



i,j=o,...,n-i n?=o(^ + i^)s;"-i 



■ A„(j b + js) 



ending the proof of equation (3.16) in the special case Zj := b + js. 

Lemma 3. With n a positive integer, a and b some complex numbers and azj + b^0, — s, 

{Zj)s;i 



det 



{azj + b)s-i 



i,j=0,...,n— 1 



{b+{n-l-j){l-a)s)^ 



A„(z). 



(3.22) 



(3.23) 



(3.24) 

-{n-2)s, 
(3.25) 
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This formula generalizes equation (C.6). 

Proof 1. From equation (2.30), with n — 1 > i > and az + b^0, —s, . . . , — (n — l)s, 

{z)s;i 1 



{az + h)s-i [az + h)s-n-i 
When a = 1, from the binomial formula (2.41) and the multiplication law (2.22) 



{z)s-i{az + b + is)a-n-i-i ■ (3.26) 



{z)s;i{z + h + is)s-n-l-i= ^ ^ j{b)s;k {z)s;i {z + is)s-n-l-i-k 

= E \ j{^)s;k{^)s;n-l-k=M^) (3.27) 

where Ili{z), a polynomial in (z)s;. of degree n — 1, is defined as in equation (3.4) with 

{0 fc = 0, 

Thus, the matrix [ci^k\i,k=o,...,n-i is triangular and its determinant is the product of its diagonal elements, 

n-l 

det[ci,fc]i,fe=o,...,n-i = W ib)s;j ■ (3.29) 
Then, when a = 1, lemma 3 follows from equation (3.3), with zj + 0, —s, . . . , — (n — 2)s, 



det 



{Zj)s;i 



{Zj + b)s;i 



0,...,n — 1 



f 1 1 

1 n <h\ r det[ni(^j)]ij=o,...,n-i 



^ (Zj + 6)s;„ 

_ib) 

(z. - 

j=0 



3 

n-l 



n^jffcrt^-W- (3-30) 



Note that for s = 0, the equation (C.6) for the power function can also be derived as above. 
When a / 1, {z)s;i {az + b + is)s;n-i-i is still a polynomial 11^(2;) of degree n— 1 in {z)s;. (or equivalently z). 
But now the evaluation of the connecting coefficients Ci^k (or even to compute det[ci^/c]i_fc=o,...,n-i we only 
need) is no longer easy since there is no simple combination law between the s-shifted factorials of z and az. 
proof 2 provides a simple proof of lemma 3 for all values of a. 

Proof 2. This proof of equation (3.25) follows the same steps as proof 2 of lemma 1. With the linear 
combinations of rows TZi — (Mi_o/Mj_i^o) "^i-i, 

M,,:=^i^ - ^;^-'^Kb+i^-m-a)s) jz.U-, , = (3.31) 

[azj + b)s-i [azo + b+{i- l)s) {azj + b) [azj + b + s)s;i-i 

Then, with azj + b ^ 0, —s, . . . , — {n — 2)s, the recurrence formula on n reads 

Ds;n{zo, • • • , Zn-i; a, b) := det [Mjj] jj^o,....n-l 



(^)(l-o)s;n-l f -pr Zj — Zq\ i , \ /o oo^ 

~1 11 >Ds-n-i{zi,...,Zn-i;a,b + s). (3.32) 



{azo + b)s-n-i V . , 
/ 1 ^ J=l 

Iteration of this equation, down to D^-i {^Zn-W a,b+ {n — l)s) = 1, ends the proof. 
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Note that in lemma 3, corresponding to a rearrangement of the factors, 



n— 1 n— 1 

n (6 + (n - 1 - - a)3)^, . = J] (6 + (n - 1 - j)s) . 

3=0 j=0 



(3.33) 



As in proof 2 of lemma 1, in the special case Zj := c + js with c some complex number, a = 1 and s 
nonzero, the determinant can be made triangular in one step, replacing Tli by the linear combination 



fe=0 



{c+{i-l)s)_ 



s:i — k 



kj {d + 2{i-l)s) 



(3.34) 



■s;i—k 



where d := b + c. Indeed, with Mij := (c + js)s:i/{d + is)s;i and d 7^ 0, — s, . . . , — (2n — 3)s, after some 
elementary algebra based on the relations (2.5)-(2.7), (2.12) and (2.13), one gets 



fe=0 



(-1) 



kJ {d+2{i-l)s)_^,._^ {d + js)s;k 

3F2(cs~^ + j, ds~'^ + -i; cs~\ ds~^ + j; 1) 



{d+{i-l)s) 

where the 3F2 is a terminating Saalschiitzian generalized hypergeometric series [12] , 

3F2(a, /?, -i; 7, 1 + a + /3 - 7 - 1) = j^f ~ ^j' i = 0, 1 



(3.35) 



Hence, 



(7)i (7 - a - /3)i 

■S* b']i - C)s;i 



{d + js),.4d+{i-l)s) 



(3.36) 



(3.37) 



which vanishes for i > j. Another proof of this identity is given in appendix D.3. Then, the determinant is 
the product of its diagonal elements and with equation (3.23), one finds 



det 



(C + js)s;» 

{d + js)s-i 



n-1 

= n 

i,j=0,...,n-l j^Q 
n-1 



= n 

3=0 



s^]\{d-c), 

{d + js)s;n- 
{d + js)s;n- 



(3.38) 



ending the proof of equation (3.25) in the special case Zj := c + js and a = 1. 



3.2. Consequences for s-shifted factorial with a negative index 

Using equation (2.24), {z)s--i = l/{z — s)-s;i, and (B.2), the following corollaries are direct consequences of 

the previous lemmas. 

Corollary 1. With n a positive integer, s some complex number and zj ^ s,2s, . . . ,{n — l)s, 



det [{zj)s-- 



J 0,...,n— 1 



det 



1 



{Zj S)_s;i 



= (-1) 



n(n-l)/2 



i,j=0,.",'"'— 1 



|n(^.)«;-(n-i)jA„(z). (3.39) 



Proof. Consequence of lemma 2. 
" See, e.g., [12] equations 2.1(30) and 4.4(3). 
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Corollary 2. With n a positive integer and s some complex number, 

1 



det 



det [{zj)s 



A„(z). 



(3.40) 



Proof. Consequence of lemma 1. 

Corollary 3. With n a positive integer, a, b and s some complex numbers and azj + b ^ s,2s, . . . ,{n— l)s. 



det 



{aZj + b)s-- 



= det 



{Zj s)-s-i 

{azj -\-b- s)-s;i 



i,j=0,...,n—l 



n 



I {b + s + {n-j){a-l)s)^,_ 



A„(z). 



(3.41) 



Proof. Consequence of lemma 3. 
Remarks: 

(i) It should be noted that equations (3.39), (3.40) and 'almost' (3.41) can be obtained, respectively, from 
equations (3.16), (3.2) and (3.25) changing for all w and i, l/{w)g-i into {w)g--i, although these quantities 
are not equal. 

(ii) A proof following the same steps as proof 2 of lemma 1 , and using the same linear combination of rows 
IZi — {Mifi/Mi--ifi)TZi-i, can also be given for corollaries 1 3. 

(iii) The extensions of lemma 1, corresponding to equations (3.3) and (3.6), apply as well to lemmas 2 and 
3 and to corollaries 1-3, see equation (C.7), e.g., with 



n-l 



fe=0 



then 



det[n,(z,)] .^^.^^^ = A det[(^j)«;-fc],.,;i=o,...,„-i 
and also, with t some complex number. 



(3.42) 



(3.43) 



det[(z,).;t-»],_^.=o,...,n-i = det[(z, +t),;_i]._^.^o_^^ 

3=0 



(3.44) 



Lemma 4. With n a positive integer and s some complex number, 

det[(.. +«;,)..-i].,.o,....„-i = ^^^^ 

Ulj=o 3-) 



A„(z) A„(w). 



(3.45) 



Proof. By the binomial formula (2.41), with Af^j := (z, + Wj)^.^-! = J2kZl Ck^) {zi)s;k (wj)s-n-i-k, the 
matrix M reads as the product of two matrices. Since the determinant of the product is the product of the 
determinants, one gets 

det[(0i + «;j)«;„_i].^.^Q = det[(^"^ ^^(zi).;fe]i,fc=o,...,n-i ^^^[('^j)^'n-i-k] j^k=o,...,n-i ■ (^-^S) 

Taking the binomial coefficients out of the first determinant and rearranging the rows of the last determinant, 
equation (3.45) follows from lemma 1. 
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Note that in all lemmas and corollaries above, the determinants considered arc anti-symmetric poly- 
nomials or rational fractions of the n variables 2:0, ... , Zn-i, therefore one expects the simplest polynomial 
alternant A„(z) to be a factor of the result. The same argument holds for A„(w) in lemma 4. 



4. Determinants with gamma functions or binomial coefficients as elements 

Using the relations (2.12)-(2.14) between the s-shifted factorial and the gamma function or the binomial 

coefficient, the results listed below arc immediate consequences of the formulae derived in section 3 with 
s = ±1. For corollaries 4-6, a direct proof following the same steps as proof 2 of lemma 1, and using the same 
linear combination of rows TZi — (Mifi/Mi-ifi) TZi-i, can also be given. It is only sketched as an example for 



corollary 4. In the special case zj 
A„(z) is given by equation (B.5). 



b + aj, with a and b some complex numbers, the product of differences 



Corollary 4. With n a positive integer, 

ri-l 

det [r{zj + i)],j^o,...,n-i = { n ^n(z) Zj ^0,-1,... 



det 



nn—i .1 
7=0 T- 



3=0 



(4.1) 
(4.2) 



Proof 1. Consequences of lemma 1. 

Proof 2. With the linear combination of rows TZi — (Mj_o/Mj_i^o) T^i-it 

Mij :=T{zj +i) {zj - zo)T{zj + i = l,...,n-l. 

Then, with Zj 7^ 0, — 1, . . ., the recurrence formula on n reads 

Dn{zQ, . . .,Zn-i) := det[Mij] 

n-l 

= r(^o) I Yl i^j ~ ^0)} -D„-l(2;i, . . . , Zn-l) . 

Iteration of this equation, down to £'i(z„_i) = T{zn-i), ends the proof of equation (4.1). 

In the special case zj := b + j ^ 0,-1, . . ., one recovers the result already published in [1] 

det [r{b + z + J)] ,,,^o,...,n-i = n + j) ■ 

j=0 



(4.3) 



(4.4) 



(4.5) 



Corollary 5. With n a positive integer, 

det ^ 



(_l)"("-l)/2 

i,j=o,...,n-i n?Jo rfe- 



A„(z) 



and for zj 7^ 0, 1, . . . , n — 2, 



det 



(?) 



i,j=0,...,n—l 



(_l)"(n-l)/2|-Q_^j^^(,) 



(4.6) 



(4.7) 



15 



See [1] equation (A. 12). 
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Proof. Consequences of lemma 2. 

In the special case Zj := 6 + j, one gets [13] 



det 



T{h + i + j) 



= (_l)"(«-l)/2 

3=0 



ij=o,...,n-i f±nb + n-l+j) 



Corollary 6. With n a positive integer and b some complex numbers, for zj ^ 0, —1, . . ., 



det 



r{z, + 1) 



T{azj + b + i) 



'^^ (6+(n-l-,)(l-a))^.r(z,-) 



and for a^;j + 6 9^ 0, 1, . . . , n — 2, 



det 



jj=0,...,n-l ''j=0 



A„(z) . 



Proof. Consequences of lemma 3. 

In the special case Zj := c + j, a = 1 and d ■.= b + c with c ^ 0, —1, . . ., one gets [13] 



det 



Tjc + i + j) 
T{d + i + j) 



;^r^-l-i 



where H^^o " c)j = U]=o {d-c + j) 
Corollary 7. With n a positive integer and zj ^ 0,-1, .. ., 

det [r(^,. - i)] . = (-l)"("-i)/2 {n r(^,- - n + 1)} A„(z) . 



Proof. Consequence of corollary 1. 
Corollary 8. With n a positive integer, 

1 



det 



i,j=0,...,n-l lVj=0^{^j) 



A„(z) . 



Proof. Consequence of corollary 2. 

Corollary 9. With n a positive integer and azj + b ^ n — l,n — 2, . . ., 



det 



T{azj +b — i) 



T{zj-i) 



n-l 

n 



T{azj + b-n+l) 



=o,...,n-i Ifio ib+'^ + in-j){a-l)) T{zj), 



A„(z), 



Proof. Consequence of corollary 3. 
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In the special case Zj := c + j and a = 1 with d := b + c n — l,n — 2, . . one gets 



det 



T{d + j-t) 
T{c + j-i) 



0,...,n— 1 



j=0 



T{d-n + l+j) 
T{c + j) 



(4.15) 



Corollary 10. With n a positive integer, 
T{zi + Wj + n - I) 



det 



r{zi + Wj) 

det 1^^ + "^ 
n — 1 



(_l)"(n-l)/2 



((n-1)!)" 



i,j=0,...,n-l (nj=oiO 
(_l)"(n-l)/2 



A„(z) A„(w) 



i,j=0,...,n-l (n"=oiO^ 



A„(z) A„(w). 



(4.16) 
(4.17) 



Proof. Consequences of lemma 4. 



5. Some examples of applications 

Let us sketch some examples of applications which motivated this work, i.e. the calculation of the probability 
density of the determinant (PDD) of random matrices. Three ensembles of n x n random matrices, with 
n = 1,2,..., have been extensively investigated, namely the orthogonal (/? = 1), unitary (/? = 2) and 
symplectic (/3 — 4) ensembles of, respectively, real symmetric, complex Hcrmitian and real quaternion self- 
dual matrices [14]. Then, the probability density of the eigenvalues x := {xj real £ Vjj = 0, . . . ,n — 1} 
reads 



Pn,/3(x) = Cn,0 |A„(x)| ]J w{Xj) 

j=0 



(5.1) 



where Cn,/3 is the normalization constant, A„(x) is defined by equation (B.l) and w{x) is a non-negative 
weight function. Quantities one computes in random matrix theory are often expressed in terms of determi- 
nants (or Pfaffians). This is the case for the expectation value of any factorized function of the eigenvalues, 
^i^) '— YYjZo Vi^j) [1'"^]- Let us show here this result only in the simplest case (3 = 2, namely with 
diJL{x) := w(x) ip{x) dx, one has 



!/x(a;o) • • • / dfi{xn- 
Jv 

^j,k-= / dii{x) Pj{x) Qk{x) 
Jt> 



i)|A„(x)| 



i! det [$ 



'•'■''^J j,fe=0,...,n-l 



(5.2) 
(5.3) 



where Pj (resp. Qk) is any monic polynomial (i.e. the coefficient of its highest power is one) of degree j 
(resp. k). Indeed, from equation (C.3), each of the two factors A„(x) can be expressed as a polynomial 
alternant and expanded as X^^g^^^o n-i} ^^P) 11^=0^ Ppj i^j)' where e{p) is the signature of the permutation 
p := {po, • • • , Pn-i}- Thereby, one gets 



n-l 

($)= Yl s{p)e{a)l[ df,{x)P,^{x)Q,^{x) 

p,o-G5„{0,...,n-l} j=0-''^ 



n-l 

n! J2 n*P.-.^- 

pG5„{0,...,ri-l} j=0 



(5.4) 



completing the proof of equation (5.2). According to the measure dfi{x) considered, one may take advantage 
of the freedom of choice of the monic polynomials in order to simplify the calculations. Thus, it may be 
useful to choose the set of orthogonal (or skew orthogonal for /3 = 1 or 4) polynomials with respect to the 
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weight [14, 15] . For example, taking $ as the identity operator, the result above with (p{x) = 1 

provides a convenient way to compute the normalization constant, e.g., for /3 = 2 



n— 1 ^ 

(C'ra,2) = n\ TT z^j '^j '■— / dxw{x)Pj{xf (5.5) 



j=0 

where Pj are the orthogonal monic polynomials for the weight w{x). 
The calculation of the PDD, 



9nAy)-= / dxo--- dx„-iP„^(3{K)6{y-xo---x„-i) 

JT) Jt> 



(5.6) 



of the random matrices we consider, is based on the use of the Mellin transform. Since this transformation 

explores a function only on the real non-negative half-axis, one needs to compute the Mellin transform of 
the restriction to y > of both the even and odd parts of the PDD, piu) ■— \ i9n fsiu) i 9n i3{—y))- From 



equations (5.6), with s some complex number, the Mellin transform of 9^ g{y) reads 



n-l 



roo 1 /" f ^ 

^nA^) ■■= / dyy'~' atfsiy) = 7; dxo--- dx^-i p„,/3(x) [] (5-7) 

Jo ^ Jv Jt> 

ip%^{x) := e^{x) \x\'-^ e+{x) := 1 e-{x) := sign{x) (5.8) 

namely, an expression of the type given by equations (5.2) and (5.3) when /3 = 2, thus 

Mi,{s) = \Cnan\ det[$±,(s)] (5.9) 
■■= dxw{x) iff Ax) Pj{x) Qk{x) . (5.10) 

Now, one can consider several ensembles of random matrices associated with the classical orthogonal 
polynomials characterized by the weight function w{x) and the domain V [14] . 

(i) For the frequently used Gaussian unitary ensemble [1] associated with the Hermite polynomials, w{x) = 
exp(— x^) with V = 71. Choosing the polynomials Pj (resp. Qk) to be the monomial x^ (resp. x'^), one finds 
[4] 18 

/■OO 
dxe-''A^{x)\xY-^ x^+'' Res>0 
-OO 

= i(i±(_i).+fe)r(£±i±^). (5.11) 

Then, the alternate elements of det[$^j,(s)]^. being zero, we can rearrange its rows and columns 

so as to collect the zero elements separate from the nonzero elements. Note that this checkerboard structure 
of the determinant is true for any w{x)ip{x) with a well-defined parity and a domain V symmetrical with 
respect to a; = 0. Thus, 

det[$+fc(s)] = det[$+._2fe(s)],-,fc=o,...,[(n-i)/2] ^^^['^%+iak+Mi,k=o,...A{n-2)/2] (5-12) 

det[$-.(^)],,,=o,...,„-r = ( (d^^t*^".-+i(^)]..^o,...,„/.)' " (5.13) 

I n odd 

See, e.g., [15] appendix A.14. 
See, e.g., [14] section 19.3. 
1^ See, e.g., [4] 6.1.1. 
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where [.t] (Icnotes the largest integer less than or equal to x. From equation (5.11), the three determinants 
above are of the type considered in coroUary 4, equation (4.5), e.g.. 



det[$+-2fe(s)],- fe^o,...,[(n-i)/2] = det[r(f + J + /c)],- ,^o,...,[(n-l 



)/2] 



[(«"l)/2] 

n j!r(i+j)- (5.14) 

j=0 



(ii) For the so-called Laguerre unitary ensemble [3], w{x) = exp(— a;) with a > —1 and I> = [0, cx)[. Still 
choosing the polynomials Pj (resp. Qk) to be the monomial (resp. x*), one finds [4] 



*ffc(s)=/ da;a;"e-^|a;|"-^a;^+'= =r(s + a + i + A;) Re s > 

Jq 

Then, again with corollary 4 equation (4.5), one obtains 



n-l 



det[*tfc(«)],,,=o,...,n-i = n^''r(^ + " + i) 

j=0 



(5.15) 



(5.16) 



the result being the same for ± since the spectrum is non-negative. 

(iii) For the so-called Gegenbauer unitary ensemble [3]: w{x) = (1 — x'^)^~^^'^ with A > i and V = [—1, 1]. 
Note that the special case A = i corresponds to the the so-called Legendre ensemble with w{x) = 1. Still 
choosing the polynomials Pj (resp. Qfe) to be the monomial x^ (resp. x''), one finds [4] ^° 



$±^(s) = j dx{l-x'^)^-^/^e^{x)\x\'-^x^+'' Res>0 



l(l±(-l)^+'=)r(A+i)^^^ 



r(^) 



s+]+k+l ^ 
2 , 



(5.17) 



Therefore, the equations (5.12) and (5.13) are still satisfied and the three determinants which occur are of 
the type considered in corollary 6 equation (4.11), e.g.. 



0,...,[(n-l)/2] 



= r(A + i) 



n[("-i)/2]-i-i 



det 



r(f +j + fc) 



r(^+A+j + A;) 



j,fc=0,...,[(n-l)/2] 



J-^'^ ,-,r(A + i+i)r(| + i) 
j- 



11 r(^+A+[(n-l)/2]+i) 



(5.18) 



(iv) For the so-called Jacobi unitary ensemble [3], w{x) = (1 — x)" (1 -l-a;)'' with a > 1, 6 > 1 and V = [—1, 1]. 
For a = b = X — ^, this ensemble is identical to the Gegenbauer ensemble above. For a ^ b, the problem is 
more complicated, in particular due to the fact that 'w{x) is no longer an even function. To illustrate the 
use of the formulae we derived, let us calculate only the normalization constant Cn,2- Choosing the monic 
polynomials Pj{x) = {x — ly and Qk{x) = (1 + x^, one finds from equation (5.3) with ^p{x) = 1 



)= j\x{ 



1 - xY+' (1 + x) 



b+k 



{-iy2 



a+b+i+j+k r(a + i+i)r(b + i + fc) 



j,k\-^ — V- V- > V - T{a + b + 2 + j + k) 

Then, the determinant in equation (5.2) is of the type considered in corollary 5 equation (4.8), 



(5.19) 



(Cn,2 



n< 



, I -Q ^_^y 2a+6+i+2j r(a + 1 + f) r(5 + 1 + j) det 

j=0 / 



1 



T{a + b + 2 + j + k) 



ri-l 



2n(n-l)-l-(a+&+l)ri j~J 



3=0 



T{a + b + n + l+j) 



j,fc=0,...,n— 1 

(5.20) 



19 See, e.g., [4] 6.1.1. 

20 See, e.g., [4] 6.2.1 and 6.2.2. 
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This result can be checked either from equation (5.5) using the constants associated with the Jacobi poly- 
nomials [16] , or from the Selberg integral [14] . 

Finally, for all these unitary ensembles (except, possibly, for the currently unknown Jacobi ensemble 
with a ^ h), the Mellin transform A4^2(''') appears to be a product, or a ratio of products, of gamma 
functions whose arguments are linear in s. Then, from the inverse Mellin transform, the PDD is expressed 
in terms of Meijer G-functions [17] . For the orthogonal and symplectic ensembles the expressions are 
more complicated [13], but we are still led to consider similar determinants. Note that, as a by-product, 
one gets also the non-negative integer moments of the PDD for g = 0, 1, . . ., 

MnAQ)--= I dygnAy)y'={^ + {-'^Y)Ml0{q + l) + {l-{-iy)M-A(l + '^)- (5.21) 
Jv 

In connection with quantum coherent states, Dr K A Penson brought our attention on the boson 
normal ordering problem, see [18-20] and references therein. Let a and a] be the boson annihilation and 
creation operators respectively, satisfying [a, af] = 1. The normal ordering of powers of boson monomials 
((at)'' a")", with n, r, s (r > s) some non-negative integers involves integer sequences of numbers which are 
generalizations of the usual Stirling numbers of the second kind, equation (2.9), and Bell numbers, whose 
values they assume for r = s = 1, 

ns ns 

{{a\y a'Y := (at)"^'--^) X^5,,,(n, k) {a\f a*^ B,,«(n) := ^ 5,,4n, k) . (5.22) 

k=s k=s 

A complete theory of these sequences of numbers has been worked out. In particular, the Br^s{n) can be 
expressed as a sum of of an infinite series of shifted factorials (generalized Dobihski formula) and moreover, 
can be considered as the n-th moments of a positive weight function Wr,s{^) with a; > 0, 

/■OO 

Br,s{n)= I dxx''Wr,s{x) ■ (5.23) 
Jo 

Extending n to complex values and using the inverse Mellin transform, one gets from above many solutions 
Wr,s{x) of the Stieltjes moment problem [19]. Generalizing this approach to the integer sequences arising 
from the normal ordering of exponentiated boson monomials, as given by equation (5.22), also provides 
solutions to Stieltjes moment problems. It happens that determinants of the type we evaluate are the 
Hankel determinants which positivity, if it can be proved, ensures the existence of the moment problem [20] . 

Let us add that the reader can find in [21] many methods of evaluations, lists of results and a wide 
bibliography on the determinant calculus. Beyond the evaluation of particular determinants, we want to 
point out that the properties of the s-shifted factorials given in section 2 emphasize similarities and connec- 
tions which exist with the power function (see another example in appendix A), thereby providing compact 
formulae and possibly a guide to find new relations. 
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Appendix A. Finite sum of s-shifted factorials of arithmetic progression 



For p a non-negative integer and a, r and s some complex numbers, we compute the finite sum of s-shifted 
factorials of arithmetic progression to n terms, 



n-l 



Zk:=a + kr Ss.p^n{a,r) ■.= '^ {zk)s;p n=l,2, ... 



(^.1) 



fc=0 

using the same trick as for the sum of powers of natural numbers. By the binomial formula (2.41), 



{A.2) 



Summing up both sides of this equation for fc = 0, . . . , n — 1 yields the recurrence formula on p, for n fixed, 

p-i 

I 



p 

;p,n(a,r-) = (^{Zn)s;p+1 - {zo)s;p+l ^ j <S's;^,n(a, Ws;p+l-«) P,n = l,2, .... (A3) 



[p + l)r 

The first two sums are independent of s. 



TliTl — 1 ) 

'S's;o,n(a,r) = n Ss;i,nia,r) = na + r. 



(^.4) 



When r = s, then Zk — s = Zk-i, and with s nonzero, an explicit expression of 5's;p_ri(a, s) can be 
obtained directly from the generalized Pascal triangle property (2.31), 



n-l 



Ss;p,n{a-, S) = ^{(Zk)s;p+1 - {Zk-l)s;p+l) 



k=0 



ip+l) 



{{Zn-l)s;p+l - {Z-l)s;p+l) 



{A.5) 



where Z-i = a — s. This result can be checked by recurrence using the general equation (A. 3). Similarly, for 
r = —s one gets 

(^.6) 



Ss;p,n{0', —s) — ^ [{zo)s;p+l — {Zn)s;p+l) ■ 



s{p + 1) 

Thus, for a = r = s = 1 one has, respectively, for the rising and the falling factorials 

{n)p+i 



S'i;p,n(l,l) = (l)p + --- + (n)p = 
5_i;p,„(l, 1) = [l]p + ■■■ + [n]p = { 



p+1 

n 



p = 



[p]p + • • • + [n]p = ^ _ T P=l,...,n 



p+1 



(^.7) 
(^.8) 







p = n+l,... . 



Further general properties follow from equations (2.5) and (2.7): 

'S's;p,n(-a, -r) = S-s;p,n{a,r) 

(X T 



(A.9) 
(^.10) 



Appendix B. Product of differences 
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With the notations of equation (3.1), the product of differences A„(z) is defined by 

rl n = l 

A„(z) := A„(2;o,...,^;„-i) := <^ JJ {zj - Zi) n = 2,3,...- 

L 0<'i<j<ri-l 

The following relations are immediately obtained with a and b some complex numbers: 
A„(6 + azo, ...,b + azn-i) = a"("-i)/2 A„(z) 

1 1 \ (-l)"("-l)/2 



A4-,---,--) = !^Ji, ^„ , A„(z) .,^0 



, ■ ■ ■ , , ""^ ) = J A„(z) b + azj^O. 

b + azo' 'b + azn-j UZoib + azj)^-^ ^' 



{B.l) 

{B.2) 
{B.S) 

(BA) 



< b + azo 

Finally, with a and b some complex numbers, in the special case Zj := b + aj, the product of differences reads 



n-l 



A„(i^6 + ai) = a"("-i)/2 JJj!. 



{B.5) 



Appendix C. Vandermonde's determinant 

ft is well known [5, 15] that the Vandermonde determinant det[(2:j)*]^^._Q is equal to the product 

of differences defined by equation (B.l), namely. 

More generally, let us consider any set of n linearly independent polynomials in z each of degree less than n. 



i{z):=J2ci,kz'' i = 0,...,n-l A := det[ci,fe]^ 7^ 0. 



{C.2) 



k=0 



Then, since the determinant of the product is the product of the determinants, one gets for the polynomial 
alternant 



det[Pi(^,)],_^.=o,...,„-i = dethfc]^ det[(z,)'=] 
= AA„(z). 



(C.3) 



Choosing the p^'s to be monic polynomials of degree i (e.g., the monomials z*), then Ci^k = for A: = 
i + 1, . . . ,n — 1 and Cj^j = 1, therefore A = 1 in equation (C.2). Now, with bi some complex numbers, it 
follows from the binomial formula that {z + biY is an other choice of monic polynomial of degree i, hence 



det[(6, + z,)^],,,=o,...,n-i=^"(^)- 



(C.4) 



When bi = b, the relation above is also a direct consequence of equations (C.l) and (B.2). 

Since (2*)"^ — {z~^y, with a and b some complex numbers, one immediately obtains from equations 
(C.l), (B.3) and (B.4) 



det 



1 



(_l)«(n-l)/2 



i,i=0,...,n-l n"Jo(^.)"-^ 



A„(z) Zjj^O 



det 



{azj + by 



^n(n-l)/2 



0,...,n— 1 



A„(z) aZj+by^O. 



(C.5) 
(C.6) 
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See, e.g., [5] section 14.311 or [15] section 7.1. 
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More generally, following the same arguments as for equation (C.3), one ean consider polynomials of the 
monomials introduced above (or even of any function), e.g., with A := det [cj^fc] ^ and azj + b 

nonzero, 

^^^^ k J nl II — i ) I''' 

Appendix D. Other proofs of equations (3.14), (3.21) and (3.35), (3.37) 

These identities can be proved by recurrence on i. Let us also give a proof which illustrates another way to 
handle shifted factorials, namely they can be generated by repeated derivations and/or integrations, e.g., 



d \3 

fv fy2 1 



1 

/o JO " - 



D.l. Other proof of equation (3.14) 

Differentiating j times (a; — 1)* x^~^^~^ in two ways (binomial formula and chain rule derivation of a product) 
[13], 



„6+j+fe-l 



fc=0 ^ ^ 

= tC){(i)V-ir}(|)'->- ,..3, 



and setting x = 1, only the term with i = i is nonzero. Thereby, one gets 

E(-l)*"' (l) [b + j + k- l]j = [b + J- {DA) 
fc=o ^ ^ 

where and thus the right-hand side, vanishes for i > j, see equation (2.10). Since from equations (2.12) 
and (2.13) 

[b + j + k- l]j = [b + i- {b + j)k {D.5) 

one recovers equation (3.14) (with s = 1 for simplicity). 

D.2. Other proof of equation (3.21) 
Assume first i > j. Then, as above, one gets 



fe=0 

22 



where the last equality is due to an overall factor y — 1 which remains after the derivation. When i < j, 
integrating j — i + 1 times {y — ly y^+^~'^ in two ways and then setting y = 1 yield 

1 



kj {b + i + k - 



py ry2 ^ 

/ dyj.i+r... dy,{yr-iyy'+'-\^^^Y.(-iy 

= / dyi (yi - 1)' yl+'-' [ dy^- - I dyj.^+i = -till b(6 + i - 1, j + 1) (£>.8) 

Jo Jyi Jyj-i U ~ 



'VI •'Vj 

where B{z,w) is the beta function [4] thus 

Now, since from equations (2.24), (2.12) and (2.13) 

(-l)-[6 + i + .-2]._,_. = ^^^^Zl),..,. (^-^^^ 

T{b + j) (-6-2(i-l))._, ib + j)k ^ • ^ 

and furthermore 

r(b + ^-i) r(6 + 2z-i) 1 

r(6 + z + j) r{b + j) {b + i-iy(b + jy ^ " ^ 

the sum over k in equation (D.7) for i > j and (D.9) for i < j does correspond to the sums considered in 
equation (3.21) (with s = 1 for simplicity). Note that since vanishes for i > j and with equation (D.IO), 
the relation (D.9) is true in all cases. 

D.3. Other proof of equations (3.35), (3.37) 

Assume first i > j. With now two variables x and y, as above, one gets 

= E(-l)'"' (fc) {c + j + k-l]j[d + i + k- 2],_,_i 

where the last equality is due to an overall factor x — 1 which remains after the derivation over x. When i < j, 
differentiating j times with respect to x and integrating j — i + 1 times over y the expression {xy— 1)* 
in two ways and then setting x = 1 and y = 1 yield 



[c + j + k-l]j 



i — k 

fe=0 



kJ {d + i + k- l)j-i+i 







/ dy2-- 


■ / dyj-i+i\ 


'yi 





x=l 



(-4^ E(-l)'"' i'l) + ~ 1]^-' B(d + i + ^ - 1, i - ^ + 1) (D.U) 



25 See, e.g., [4] 6.2.1 and 6.2.2. 
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where B(z,?ii) is the beta function. Thereby, after some elementary algebra, using equations (2.24), (2.12) 
and (2.13), one gets 



i — k 



i\ [c + j + k- l]j- 
kj (d + i + fc - l)j_j+i 




{d + i - l)i[c + i - 



r(d + 2z-i) r(c + j) mid~c). 



{D.15) 



r{d + j) r{c + i) id + j),{d + i-i). 



where the last equality follows from the binomial formula (2.46). Now, since 




{D.16) 



r{d + 2i-l) T{c + j) [c + i-l],-k (c + j)fc 



{D.17) 



T{d + j) r(c + z) [d + 2{i-l)],^k {d + j)k 



the sum over k in equation (D.13) for i > j and (D.15) for i < j does correspond to the sum over k in 
equation (3.35) (with s = 1 for simplicity). Note that since vanishes for i > j and with equation (D.16), 
the relation (D.15) is true in all cases. 
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